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The  process  o£  wiping  a  thickened  chemical  contaminant  from 
a  surface,  where  it  has  been  inadvertantly  transferred  and  may 
pose  a  safety  hazard,  is  investigated  with  regard  to  the 
fundamental  fluid  mechanics.  Of  special  interest  is  the  effect 
viscoelastic  fluid  behavior  has  on  the  process  and  specifically 
the  amount  of  fluid  remaining  after  wiping  .  Rigid  aluminum 
blades  of  variable  geometry,  accurately  positioned  above  a  smooth 
rotating  cylindrical  surface,  provide  an  experimental  model  of 
the  process.  Residual  fluid  thickness  data  as  a  function  of 
blade  height  for  newtonian  liquids  are  compared  to  a  simple 
lubrication  theory  model,  while  a  more  thorough  Finite  Element 
solution  to  the  problem  is  planned  for  the  future.  Preliminary 
results  indicate  the  potential  to  separate  shear  thinning  and 
elastic  effects  experimentally  through  the  proper  choice  of  test 
fluids. 


The  process  of  wiping  a  contaminating  fluid  from  a  surface 
poses  an  interesting  fluid  dynamics  problem  and  is  of  practical 
interest  as  well.  Although  a  wiping  step  alone  is  rarely 
sufficient  to  render  a  surface  "clean"/  optimizing  the  wiping 
process  (minimizing  the  amount  of  contaminant  remaining)  can 
considerably  reduce  the  requirements  of  additional/  more  costly/ 
cleaning  processes, 

A  simple  geometry  for  investigating  the  wiping  process  is 
shown  in  Figure  1.  Hsu  (1984)  considers  a  similar  geometry  to 
study  the  the  forces  generated  in  a  blade  coating  process.  His 
analysis  is  the  basis  for  a  simple  analytical  model  presented 
here,  with  two  modifications.  First,  his  primary  interest  is  in 
the  pressure  field  generated  under  the  blade,  while  we  are 
concerned  with  the  amount  of  liquid  left  behind  after  the  wiping 
process.  Secondly,  he  considers  only  a  specific  case  of  our  more 
generalized  geometry,  leaving  open  the  question  of  how  sensitive 
the  results  are  to  blade  orientation  and  placement.  Following 
Hsu,  a  perturbation  method  is  used  to  obtain  an  analytical 
solution  for  newtonian  liquids,  and  justification  is  given  for 
retaining  only  the  zeroth  order  solution  (equivalent  to 
lubrication  theory  -  Cameron,  1966).  A  specially  designed  wiping 
apparatus  is  used  to  collect  residual  fluid  thickness  data,  for 
comparison.  Data  for  newtonian  liquids  show  good  agreement 

with  the  model  under  some  conditions,  and  provide  insight  into 
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the  limitations  of  the  model  under  others. 

The  model  is  extended  to  purely  viscous  non-newtonian 
liquids,  by  incorporating  a  shear  dependent  viscosity  function 
into  the  equations.  Experiments  with  viscoelastic  liquids 
demonstrate  that  the  residual  fluid  thickness  can  differ  markedly 
from  values  predicted  by  lubrication  theory  for  purely  viscous 
fluids . 
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Newtonian  Liquids? 


Figure  2  gives  a  definition  sketch  for  the  geometry  and 
coordinates.  We  start  with  the  steady  state  Navier-Stokes 
equations  (neglecting  gravity)  in  two  dimensions,  written  in 
cartesian  coordinates  as 


(1) 


The  continuity  equation  is 
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We  nondimensional ize  the  variables  with  a  set  of  scaling 
factors  as  follows; 
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the  specific  system  (the  wiper)  of  interest,  but  specify  that 
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We  define  a  Reynolds  number  (Re)  as 

Re  =  P-iLb' 

>«-L* 


and  a  geometrical  parameter  o(  as 

'  Lx 

Then  Eqs.  1-3  become 
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It  is  useful  to  rewrite  Eg.  6  as 
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Boundary  conditions  for  the  wiper  take  the  forms  of  the  following 
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equations  (see  Fig,  2) 
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Utilizing  the  boundary  conditions  on  H5  we  may  replace  Eq.  7 
with  an  equivalent  continuity  expression  o£  the  £orm 

4>  -  X  (12) 

v(t> 

It  is  not  difficult  to  see  that  the  constant  X  defined  above  is 
a  dimensionless  residual  fluid  thickness.  The  proportionality 
depends  upon  the  choice  of  the  length  scale  £  used  in  defining 

1  • 

Also  of  interest  is  the  separating  force  (per  unit  width 
transverse  to  the  flow) ,  defined  nondimensionally  as 


The  goal,  then,  is  to  solve  Eqs.4,  5,  12,  and  13  for  \  and  §"  as 
functions  of  o<  ,  Re  ,  and  the  shape  factors  that  enter  into  the 


description  of  the  nip  region  through  the  specific  forms  of  CT(^^ 
and  JSU'}  . 

So  far,  the  lubrication  approximations  have  not  been 
imposed.  Equations  8  and  9  express  the  no-slip  assumption  on  the 
kinematics.  Boundary  conditions  on  pressure  have  been  the 
subject  of  numerous  discussions  which  will  not  be  repeated  here. 

We  select  the  simplest  choice  (the  Swif t-Steiber  conditions 
-  Dowson  and  Taylor,  1979),  in  Eq.  10  and  11,  and  we  note  that 
the  positions  and  correspond  to  the  ends  of  the  blade  (XI 
and  X2).  By  its  absence  from  the  boundary  conditions  we  see  that 
surface  tension  effects  are  assumed  to  be  absent. 

An  analytical  solution  to  Eqs.  4-6  may  be  found  by  a 
perturbation  method. The  perturbation  parameters  are  o< 
(accounting  for  curvature  effects)  and  Re  (for  inertial  effects). 
Each  of  the  variables  is  expanded  as  follows: 
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Variables  with  superscript  (00)  give  the  zeroth-order  lubrication 
approximation.  Those  with  (10)  and  (01)  give  the  first  order 
corrections  for  curvature  and  inertia,  respectively. 
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( Lubrication)  Solution 
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When  Eqs.  14  are  substituted  into  Eqs.  4,  5,  and  12,  and 
terms  in  and  Re  are  dropped,  we  recover  the  lubrication 
equations  in  the  form 
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Boundary  conditions  8-11  hold,  but  do  not  include  any 
constraint  on  ^  ,  which  does  not  appear  in  this  level  of 
approximation. 

The  solutions  may  be  written  as  follows: 
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The  results  are  general  (to  this  level  of  approximation)  and 
can  be  made  particular  upon  a  choice  of  geometry,  and  O'(T)  , 
and  kinematic  boundary  condition  c|>(^ • 

Results  presented  by  Hsu  indicate  that  the  higher  order 
corrections  may  not  be  necessary  for  our  purposes.  First  order 
corrections  for  curvature  contributed  less  that  0.5%  to  the 
residual  fluid  thickness  (X),  in  the  specific  case  he  examined, 
and  the  zeroth-order  solution  was  within  1%  of  a  numerically 
generated  solution  of  the  Navier-Stokes  equations  (Eqs.  1  and  2). 
First  order  effects  would,  however,  be  important  if  our  primary 
concernwas  for  information  on  blade  loading  ($•).  Hsu's  analysis 
demonstrated  that  first  order  corrections  contributed  as  much  as 
10%  to  the  dimensionless  blade  loading.  Without  attempting  to 
derive  the  higher  order  solutions,  we  examine  the  results 
obtained  from  the  zeroth-order  approximation,  without  relying  too 
heavily  on  the  quantitative  predictions  of  blade  loading. 

The  particular  choice  of  geometry  and  kinematic  boundary 
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condition,  as  shown  in  Figure  2,  are 


1 i  ' 


Roll  Surface  —  —  -  (23a) 

Blade  Surface  ^  —  V\  -  x""^av\  (tS)  (24a) 

where  H  *  ^(0)  -  Ha  -  X0  Tan(  u?  ) 

H  is  the  smallest  separation  (in  the  y-direction) 
between  the  blade  and  roll. 

i.e.  H„  -  min{  ^'(x)  -  <j'(x)}  -  ^'(XO)-  <l'(X0) 

X0  »  location  of  minimum  blade  separation 
i.e.  X0  ■  min{  X2,  max{  XI,  R  sin(  to  )  }  } 


X-Velocity 


Y-Velocity  v\  -  V0  -  -  U  ( 


In  dimensionless  form. 


(25a) 


(26a) 


v\  T o.'va  ( 

£  ‘  S 


(24b) 


(25b) 


(26b) 


An  obvious  choice  foe  the  length  scales  L v  and  h  ,  is 

L*  “  R 
h  -  H 


where  H  is  the  blade  height  at  x-0.  However,  previous  results 
are  presented  with  LK  «  (r  h)land  for  the  purpose  of  comparison 
is  used  here  as  well.  The  definition  of  h  requires  a  closer 
look.  In  addition  to  h  ,  three  parameters  are  required  to 
specify  the  geometry?  the  blade  length,  L,  blade  angle,  ,  and 
the  blade  displacement, X2.  Hsu  considered  the  case  of  w»  o  and 
X2  ■  0,  for  which 

pw*  c 


and 

H  -  H„ 

where  the  obvious  choice  for  h  is  H# .  In  this  case,  H„ 
corresponds  to  the  blade  height  at  the  downstream  edge  of  the 
blade  (X 2)  and  the  minimum  separation  between  the  roll  and  blade. 


11 


When  X2  »  0  an  appropriate  choice  of  reference  length  ,  h,  is 
not  obvious.  Three  possibilities  are: 

1)  Separation  between  the  blade  and  roll  at  the 
downstream  edge  of  the  blade. 

h  -  ( X 2 )  -  G*  (JI2) 

2)  Separation  at  X  -  0. 

h  -  H 

3)  Minimum  separation  between  blade  and  roll 

h  -  H 

o 

For  this  study,  the  preferable  reference  length,  h  ,  is  the 
blade  height  at  X  ■  0,  for  X2  ^  0,  and  at  the  downstream  edge  of 
the  blade  for  X2  <  0.  The  downstream  edge  is  chosen  when  X2  <  0 
to  avoid  negative  values  of  h  .  Otherwise,  the  choice  of  h  ■  H, 
for  X2  £  0  ,  simplifies  Eq.  24a  nicely.  The  final  expressions 
defining  the  geometry  and  kinematic  boundary  conditions  are: 

For  X2  }  0  : 

<?(%')  -  'it1  (230) 

(S  (  -  I  -  T*"  (24c) 


(25c) 


results 


figures  3  and  4  contain  justification  for  considering  the 
zeroth  order  perturbation  solution  without  higher  order 
contributions  for  curvature.  Hsu's  results  for  X  as  a  function 
of  H#/R  indicate  that  first  order  corrections  are  small  relative 
to  the  sensitivity  of  \  to  variations  in  blade  angle  and 
displacement.  In  the  geometry  chosen  by  Hsu,  for  example,  a  two 
degree  error  in  blade  angle  or  2  mm  displacement  of  the  19  mm 
blade,  results  in  up  to  10%  change  in  X  compared  to  less  than 
0.5%  for  first  order  curvature  effects.  This  is  not  to  say  that 
first  order  effects  are  always  insignificant.  When  blade  loading 
is  of  primary  interest,  as  it  is  for  Hsu,  they  are  important  and 
should  be  included.  This  is  shown  in  Figures  5  and  6  where 
variations  in  blade  loading  associated  with  curvature  effects  are 
comparable  to  variations  due  to  changes  in  blade  angle  and 
displacement. 

The  results  are  very  sensitive  to  blade  angle  as  shown  in 
Figures  3  and  5.  It  is  unlikely,  in  practice,  that  the  blade 
angle  can  be  set  to  within  one  degree.  It  is  interesting  to  note 
that  in  the  limit  of  zero  blade  height,  the  limiting  value  of  X 
is  2/3,  for  a  zero  degree  angle  and  becomes  one  for  positive 
blade  angles,  regardless  of  how  small.  Increasing  the  blade 
angle  increases  the  residual  fluid  thickness  and  decreases  the 
blade  loading  for  a  given  blade. 


Figures  4  and  6  show  the  effect  of  blade  displacement.  The 
use  of  blade  displacement  as  a  parameter  is  somewhat  artificial 
because  by  re-defining  the  reference  frame,  displacements  can  be 
accounted  for  by  appropriate  changes  in  blade  angle.  By  defining 
the  y-axis  as  the  radius  that  passes  through  the  downstream  edge 
of  the  blade,  the  displacement  can  be  eliminated  as  a  parameter. 
The  reason  displacement  is  considered  in  this  analysis  is  for  the 
the  benefit  of  providing  a  sensitivity  analysis  for  an 
experimental  system  with  a  fixed  coordinate  system,  like  that 
used  by  Hsu  and  described  later  in  this  report.  Notice  the 
asymmetry  apparent  in  Figure  6.  Displacing  the  blade  upstream 
(negative  x-direction)  causes  a  significant  decrease  in  blade 
loading,  while  positive  displacements  cause  less  of  an  effect  in 
the  opposite  direction.  The  blade  loading  has  a  maximum  at 
X2-L/2  for  a  zero  degree  angle,  corresponding  to  the  geometry  for 
which  the  point  of  tangency  of  the  blade  is  located  at  the  center 
of  the  blade.  Positive  displacements  cause  the  blade  loading  to 
increase  until  X2-L/2,  after  which  further  positive  displacement 
will  decrease  the  blade  loading. 

Figures  7  and  8  show  the  effect  of  blade  length  on  \  and  . 
The  range  of  blade  length  over  which  the  lubrication  theory 
results  are  meaningful  is  limited  by  two  major  assumptions. 
Lubrication  theory  requires  that  entrance  and  exit  effects  are 
negligible,  L/H>>1,  arid  that  the  blade  and  roll  be  nearly 
parallel  to  minimize  curvature  effects,  L/R<<1.  When  the  blade 
is  large,  relative  to  the  roll  radius,  it  extends  so  far  that  the 
curvature  of  the  roll  is  significant  and  the  parallel  flow 
assumption  is  violated.  At  the  other  extreme,  solutions  for 


15 


small  blades  are  only  valid  for  very  small  blade  heights,  because 
as  the  blade  height  increases,  the  entrance  length  also  increases 
and  entrance  effects  become  significant. 


Non-Newtonian  Liquids: 


Many  of  the  liquids  of  interest  are  polymeric  solutions  or 
melts  which  exhibit  non-newtonian  behavior.  We  anticipate  that 
the  response  of  such  solutions  to  the  wiping  process  is 
complicated  by  non-linear  elastic  and  viscous  effects,  which  can 
be  extremely  difficult  to  model.  As  a  first  step  in 
understanding  these  phenomena,  we  turn  to  an  analysis  of  the  role 
of  non-linear  viscous  effects  by  incorporating  a  simple  shear 
dependent  viscosity  function  into  the  analysis. 

We  begin  with  the  steady  two-dimensional  dynamic  equations 
in  the  form 


(26) 


(27) 
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along  with  the  continuity  equation 


Vu 

t>x 


(28) 


A  useful  model  foe  shear  viscosity  is  the  power  law#  a 
purely  viscous  (i.e.  non-elastic)  constitutive  equation  that  we 
write  as 

x  *  u  A  <29> 


where  the  rate  of  deformation  tensor  has  two-dimensional 
cartesian  components 


The  viscosity  function  is  written  as 

/U  -•  K  ( i  \'y 

where 


(30) 


(31) 


(32) 


As  in  the  newtonian  case,  the  variables  are 
nondimensional ized  and  the  dynamic  equations  display  terms  to 
zero  and  first  order  in  *<  and  Re.  The  variables  are  then 
written  as  in  Eq.  14,  and  the  solutions  to  various  ordersin 
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and  Re  may  be  obtained. 

The  nondimensionalization  is  identical  to  that  used  in  the 
newtonian  case ,  except  that  the  newtonian  viscosity  constant 
of  those  equations  is  everywhere  replaced  by  a  factor 


v  I  U  \  *’v 

M  =  (33) 

We  refer  to  M  as  the  apparent  viscosity  at  the  (nominal)  shear 
rate  U/H0. 

In  dimensionless  form,  Eq.  31  becomes 

-»)'}* 

and  itself  is  included  in  the  perturbation  expansion  in  the 
form 


As  in  the  newtonian  case,  we  concern  ourselves  only  with  the 
zeroth  order  solution,  but  for  different  reasons.  Although  the 
power  law  constitutive  equation  is  a  useful  mathematical  model 
for  investigating  the  effect  of  a  shear  dependent  viscosity,  it 
is  not  expected  to  model  quantitatively  the  behavior  of  any 
"real"  fluid  involved  in  the  wiping  process.  An  immense  amount 
of  algebra  is  required  to  generate  the  first  order  contributions, 
and  it  is  more  reasonable  to  extract  information  from  the  zeroth 
order  solution  and  compare  that  to  the  behavior  of  some  real 
fluids. 

The  Zeroth  Order  Solution: 


The  equations  are  (Cf:  Eqs.  15-17) 
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where 
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Boundary  conditions  8-11  hold. 

Unfortunately,  no  analytical  solution  of  these  equations  may 
be  found.  Instead,  we  adopt  a  procedure  due  to  Horowitz  and 
Steidler  (1960)  and  produce  numerical  solutions  for  the  variables 
of  interest. 

One  advantage  of  the  Horowitz/Steidler  procedure  is  that  it 
can  be  generalized  to  any  purely  viscous  constitutive  equation. 
We  rewrite  Eq.  36  as 


where 
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but  jj.**  is  left  unspecified  at  this  stage, 
function  £(£***)  as 

r. 


We  define  a 


(42) 
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(in  the  following,  we  drop  all  superscripts  that  denote  the 
zeroth-order  solution.) 

Integration  of  Eq.  40  yields 


where 
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Integration  of  Eq.  42  using  Eq.  9  gives 


<P  =  J  jl^  Af  + 


Using  Eq.  43,  we  write  this  as 
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With  the  boundary  condition  cf>  *  0  at  mb(t')  we  may  write  this 
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where  we  define 
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With  these  results  we  may  write  \  in  the  form 
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where 
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(49) 


From  Eq.  43  we  write 
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If  we  regard  ,  X  ,  and  "TT  as  the  unknowns,  we  may 

use  Eqs .  46,  48,  and  50  in  an  iterative  approach  as  follows: 


A  constitutive  equation  is  selected/  so  that  £(x) is 
given  as  an  explicit  function.  Hence  the  integrals 
F  and  G  may  be  obtained. 

A  trial  value  for  X  is  assumed. 

Regarding  ,  T^.  ,  and  TT  as  the  unknowns  at  each 
position  |  ,  we  assume  a  t  value,  calculate  cr(5;), 

,  and  and  solve  for  ,  and  IT  at  each 

5  on  the  range  £  <  $  <  ?x  • 

i 

From  TT  and  the  boundary  conditioner  »  0  at  5  ■?, 
we  integrate  to  obtain  the  function  TV  (5).  If 
does  not  satisfy  the  boundary  condition  Tr ■  o  at 
f  ,  we  change  \  and  carry  out  another  cal¬ 

culation  of  ,  until  satisfactory  convergence 
is  attained. 


This  method  is  checked  by  carrying  out  the  procedure  for  a 
newtonian  fluid,  for  which  an  analytical  solution  for  TTtf)  and 
\  is  available  (through  the  zeroth-order  newtonian  analysis 
presented  earlier).  Satisfactory  results  are  obtained,  as  shown 
in  Figures  9a  and  9b, 
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RESULTS 


For  a  power  law  fluid,  for  which 
-!  -  \ 
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one  obtains  the  results  shown  in  Figures  10-17  (See  Table  2  for 
Key  to  Figures  10-17).  <£  is  calculated  as  before  (i.e. using 
Eq,  13),  but  for  a  power  law  fluid 


$  =  (  li)[  nuR  ^  154) 


with  M,  the  apparent  viscosity,  calculated  from  Eq.  33.  We  note 
that  U/H  Q  is  an  arbitrary  but  useful  definition  of  a 
characteristic  shear  rate  in  this  geometry. 

Figures  10-12  and  14-16  show  that  the  residual  fluid 
thickness  increases  as  the  liquid  becomes  increasingly  shear 


thinning  (i.e.  as  n  decreases)  and  is  relatively  unaffected  by 


small  changes  in  blade  angle  or  displacement.  The  increase  is 
small  and  shows  a  slight  dependence  on  H/R. 

Figures  13  and  17  show  the  effect  of  shear  thinning  on  blade 
loading  with  blade  angle  and  displacement  as  parameters, 
respectively.  Blade  height  is  predicted  to  have  a  significant 
effect,  but  small  changes  in  angle  and  displacement  do  not. 
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EXPERIMENTAL 


Apparatus : 

Figure  18  shows  a  sketch  of  the  wiping  apparatus.  A 
cylindrical  steel  roll  is  half  immersed  in  the  test  fluid  and 
driven  at  constant  speed/  U.  Liquid  is  entrained  by  the  moving 
roll  and  is  partially  wiped  by  the  blade/  while  the  remainder 
emerges  downstream  as  residual  fluid  of  thickness/  h.  The  blade 
is  held  in  a  rig  which  permits  accurate  positioning  of  the  blade 
relative  to  the  roll. 

The  roll  has  a  radius  of  2.75  inches/  accurately  machined  to 
within  0.5  mil/  and  an  axial  length  of  5  inches.  Blades  made 
from  aluminum  are  interchangeable  and  have  the  geometry  defined 
in  Figure  19.  Table  2  summarizes  the  specific  geometries 
discussed  in  this  report. 

Residual  fluid  thickness  is  measured  by  a  direct  contact 
method  (Greener/  1978)  using  a  micrometer  driven  needle 
positioned  60  degrees  downstream  of  the  blade.  The  contact  of 
the  needle  with  the  surface  is  observed  through  a  microscope. 


Procedure: 

It  is  unfortunate  that  in  the  limit  of  very  small  blade 
height  (Ho-»0),  when  the  lubrication  approximations  are  most 
likely  to  be  valid,  good  experimental  data  are  very  difficult  to 
obtain.  Because  of  unavoidable  experimental  uncertainty  and  the 
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finite  precision  with  which  the  roll  is  machined,  there  is  always 
some  error  associated  with  measuring  both  the  blade  height,  H, 
and  residual  fluid  thickness,  h.  The  small  but  finite 
uncertainties  amount  to  large  relative  errors  in  residual  fluid 
thickness  as  0.  If  the  blade  height  and  the  corresponding 
residual  fluid  thickness  are  accurately  known  at  one  setting, 
however,  precise  incremental  changes  can  be  made.  A  procedure 
for  gathering  and  presenting  data  is  adopted,  which  minimizes  the 
experimental  uncertainty  by  taking  advantage  of  the  precision 
involved  in  measuring  increments  in  H.  and  h. 

Blade  height  is  arbitrarily  set  to  a  large  value  (greater 
than  50  mil)  and  the  corresponding  residual  fluid  thickness  is 
measured.  The  blade  height  is  then  incrementally  decreased,  with 
the  associated  micrometer  reading  for  residual  fluid  thickness 
recorded  for  each  increment,  until  the  blade  touches  the  roll 
(corresponding  to  H% ■  0).  The  roll  is  then  stopped  and  the 
needle  is  lowered  until  it  too  touches  the  roll,  and  the 
micrometer  reading  is  recorded  (corresponding  to  h  -  0).  The  data 
generate  a  curve  of  h  as  a  function  of  H  .  This  curve  as  a  whole 
can  then  be  shifted,  within  bounds  set  by  the  experimental 
uncertainty  (a  maximum  of  1  mil  in  any  direction),  so  that  it 
passes  smoothly  through  the  origin.  In  this  way,  the  errors 
involved  in  individual  data  points  are  averaged  out  over  the 
entire  curve. 


RESULTS 


The  range  over  which  certain  experimental  variables  are 
considered  is  determined  by  experimental  limitations.  In  order 
to  flood  the  inlet  (upstream)  region  of  the  blade,  a  sufficient 
quantity  of  fluid  must  be  entrained  from  the  reservoir  as  shown 
in  Figure  2a.  The  amount  entrained  depends  upon  both  the  roll 
speed  and  the  fluid  properties  (Campanella,  1984).  For  a  given 
fluid,  the  roll  speed  is  limited  on  the  low  end  by  this  flooding 
criterion.  On  the  high  end,  there  is  also  a  limit,  corresponding 
to  the  onset  of  instabilities  and/or  air  entrainment  ,  Under 
unstable  wiping  conditions,  the  residual  fluid  layer  is  no  longer 
smooth  and  often  develops  a  "ribbed"  pattern  (Pearson,  1960, 
Sullivan,  1979)  transverse  to  the  direction  of  flow.  Positive 
blade  placements  and  negative  blade  angles  are  found  to  promote 
this  instability.  When  stability  is  not  a  problem,  there  is  a 
limit  on  roll  speed  above  which  air  is  entrained  into  the  test 
fluid  at  the  interface  where  the  roll  surface  re-enters  the  bath 
(Burley,  1984).  The  range  of  roll  speeds  investigated  is  from  2 
to  20  in/sec,  depending  on  the  test  fluid.  To  obtain  flooding  of 
the  blade,  a  lower  limit  on  the  viscosity  of  newtonian  test 
fluids  is  about  3  poise. 

Two  newtonian  liquids  are  examined:  glycerine  and  Karo 
syrup.  Under  conditions  that  the  blade  inlet  is  flooded,  no 
effect  of  linear  roll  speed  on  residual  fluid  thickness  is 
observed.  Figure  20  shows  data  for  glycerine  and  Karo  syrup  with 
a  1/4  in.  blade  at  0  degrees  and  the  corresponding  results 
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predicted  by  lubrication  theory.  Results  for  glycerine  and  Karo 
syrup  are  identical  (within  experimental  uncertainty)  even  though 
their  viscosities  differ  by  almost  an  order  of  magnitude.  In 
other  geometries,  similar  agreement  between  the  two  is  found. 
The  remainder  of  the  newtonian  results  presented  are  for 
glycerine,  with  the  assumption  that  they  are  representative  of 
newtonian  liquids  in  general. 

In  Figure  20,  agreement  between  the  data  and  lubrication 
theory  is  good  for  small  blade  heights  ,  but  becomes  poor  as  H, 
increases,  while  in  Figure  21,  the  data  and  theory  are  in  good 
agreement  up  to  the  maximum  blade  height  shown  (50  mils).  The 
poor  agreement  with  small  blades  at  large  blade  heights  is  due  to 
entrance  effects  not  accounted  for  by  lubrication  theory.  When  Ho 
■  50  mils,  L/H  «  10  for  the  1/2  in.  blade,  but  only  5  for  the  1/4 

O 

in.  blade.  Data  indicate  that  L/H#  >  10  is  a  reasonable  estimate 
of  the  requirement  that  entrance  effects  are  negligible.  Notice 
that  when  entrance  effects  are  significant,  data  lie  above  the 
theory  as  expected  on  the  basis  of  momentum  transfer  arguments 
(Appendix  A).  Blades  larger  than  1/2  in.  are  not  used  for  two 
reasons.  With  the  1/2  in.  blade  and  2.75  in.  radius  roll,  L/R  » 
0.18.  Larger  blades  would  likely  invalidate  the  negligible 
curvature  assumption,  L/R  <<  1,  and  at  the  same  time  require  a 
large  amount  of  fluid  to  flood  the  inlet  region  of  the  blade. 
This  flooding  problem  (see  Figure  2a),  limits  the  blade  geometry 
in  a  number  of  cases. 

Changing  the  blade  angle  from  zero  to  10  degrees  affects  the 
residual  fluid  thickness  as  shown  in  Figures  22,  23,  and  24,  for 
three  blade  lengths.  Surprisingly  good  agreement  between 
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experimental  data  and  simple  lubrication  theory  is  found  for  the 
1/2  in.  blade.  Note  that  as  the  blade  size  decreases,  the  data 
fall  further  above  lubrication  theory  predictions  due  to  entrance 
effects.  Blade  angles  larger  than  10  degrees  are  not  discussed 
because  of  large  curvature  effects  and  the  problem  of  flooding 
the  inlet  region  of  the  blade.  Some  experiments  performed  with 
large  blade  angles  (e.g.  45  degrees)  produced  data  that  fell 
below  lubrication  theory  predictions,  but  it  is  unknown  whether 
this  effect  was  due  to  curvature  or  an  inadequately  flooded  inlet 
region. 

Figures  25  and  26  show  the  effect  of  blade  displacement  and 
are  also  good  examples  of  the  above  mentioned  flooding  problem. 
Once  again,  data  for  the  largest  blade  <L«l/2  in)  are  in  good 
agreement  with  theory,  but  only  for  three  of  the  four 
displacements  shown.  The  data  for  the  largest  negative 
displacement  (X2--0.5)  fall  below  the  theory.  This  is  a  prime 
example  of  not  enough  fluid  being  entrained  by  the  roll  to  flood 
the  inlet  region  of  the  blade.  The  problem  is  exacerbated  when 
the  blade  angle  is  increased,  as  shown  in  Figure  26. 


Non-Newtonian  Liquids: 

Predicting  the  behavior  of  viscoelastic  liquids  is  difficult 
in  even  the  simplest  flow  problems  and  is  complicated  further  by 
the  converging-diverging  nature  of  the  wiping  process.  We  do  not 
expect  good  quantitative  agreement  between  the  simple  power  law 
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lubrication  theory  model  and  residual  fluid  thickness  data  for 
real  viscoelastic  fluids,  but  rather  hope  to  gain  insight  into 
the  effect  the  shear  dependent  viscosity  has  on  the  results. 

Figure  27  shows  viscometric  data,  taken  with  a  Rheometrics 
Fluids  Rheometer,  for  a  shear  thinning  carboxymethy  lcel lulose 
(CMC)  solution.  For  a  linear  roll  speed  of  10  in/sec  and  blade 
heights  ranging  from  10  to  50  mils,  the  nominal  shear  rate  ranges 
from  200  to  1000  1/sec.  In  this  range  the  power  law  index,  n,  is 
about  0.6.  Note,  however,  that  the  fluid  also  exhibits  a  finite 
normal  stress.  At  a  representative  shear  rate  of  500  1/sec,  the 
recoverable  shear,  ,  defined  by 
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is  also  about  0.6.  The  recoverable  shear  is  a  parameter  commonly 
used  to  estimate  the  relative  importance  of  elastic  and  viscous 
effects.  When  we  examine  data  for  the  CMC  solution,  we  observe 
the  results  of  both  shear  thinning  and  elastic  effects,  while  the 
power  law  lubrication  theory  only  accounts  for  the  shear  thinning 
behavior . 

Figure  28  shows  that  the  residual  fluid  thickness  of  the  CMC 
solution  is  considerably  greater  than  that  for  glycerine  and  is 
also  greater  than  predicted  by  lubrication  theory.  Since  shear 
thinning  behavior  is  accounted  for  in  the  model,  it  is  reasonable 
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to  suspect  that  the  increase  in  thickness  is  due  to  elastic 
effects.  To  test  this  hypothesis,  small  amounts  of  high 
molecular  weight  polyacrylamide  (PAM)  are  added  to  the  CMC 
solution,  to  increase  the  recoverable  shear  without  significantly 
affecting  the  power  law  index.  Data  for  two  such  solutions  are 
included  in  Figure  28,  one  with  a  recoverable  shear  of  0.9  and 
the  other  with  1.4.  Note  that  by  increasing  the  recoverable 
shear,  the  residual  fluid  thickness  decreases  rather  than 
increases.  The  data  for  the  solution  with  the  largest  recoverable 
shear  fall  below  those  of  the  newtonian  fluid.  This  iB  opposite 
to  the  effect  predicted  by  the  power  law  lubrication  theory 
model. 

In  an  attempt  to  isolate  elastic  and  shear  thinning  effects, 
the  behavior  of  an  elastic  fluid  with  insignificant  shear 
thinning  (i.e.,  a  Boger  fluid  -  Boger,  1977)  is  examined.  The 
fluid  is  made  by  adding  a  small  amount  of  high  molecular  weight 
PAM  to  corn  Byrup.  Figures  29  and  30  show  viscometric  data  for 
the  Boger  fluid  as  well  as  for  a  CMC  solution  and  glycerine. 
Note  that  the  viscosity  of  the  Boger  fluid  is  nearly  constant 
over  the  shear  rate  range  for  which  data  were  collected.  The 
Boger  Fluid  also  exhibits  large  normal  stresses  in  steady  shear 
relative  to  the  highly  shear  thinning  (n«0.6)  CMC  solution. 

Figure  31  shows  a  drastic  decrease  in  residual  fluid 
thickness  with  the  Boger  fluid  relative  to  glycerine.  Note  that 
although  the  CMC  data  fall  above  the  power  law  lubrication  theory 
curve,  it  is  in  better  agreement  for  the  larger  1/2  in.  blade 
than  for  either  the  1/4  in.  blade  (Figure  32)  or  the  0.068  in. 
blade  (Figure  33).  This  is  attributed  to  entrance  effects,  which 
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for  flow  in  pipes  are  known  to  be  considerably  linger  for 
viscoelastic  fluids  than  those  encountered  with  newtonian  fluids 
(Tung,  1978).  Results  in  Figures  31-33  demonstrate  that  the 
elastic  solution  exhibits  a  greatly  decreased  residual  fluid 
thickness  in  comparison  to  lubrication  theory  models.  The  amount 
of  viscoelastic  fluid  left  behind  after  the  wiping  process 
appears  to  depend  on  the  competing  effects  of  shear  thinning, 
causing  it  to  increase,  and  elasticity,  causing  it  to  decrease. 
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Experimental  results  for  newtonian  liquids  are  very  useful 
for  determining  the  range  of  validity  of  the  lubrication  theory 
model,  but  unfortunately,  that  range  is  severely  limited.  One 
obvious  inadequacy  of  the  model  is  the  negligible  curvature 
assumption,  which  is  only  valid  for  very  small  wiper  angles. 
Another  is  the  assumption  of  negligible  entrance  effects,  which 
appears  to  be  violated  for  thin  wipers.  The  model  also  fails  to 
take  into  account  surface  tension  effects  by  using  the  Swift- 
Stelber  boundary  conditions.  A  useful  estimate  of  the  magnitude 
of  surface  tension  effects  is  the  dimensionless  capillary  number 
defined  by 


K)  ,  AX 


The  fact  that  N  is  of  order  unity  in  some  of  our  experiments 

Co. 

indicates  that  surface  tension  may  be  significant.  Gravitational 
effects  are  neglected  in  the  model  as  well,  and  can  be  estimated 
by  another  dimensionless  group,  the  Stokes  parameter,  defined  by 


St  •  Ui* 


If  the  blade  height  is  large  and  the  roll  speed  is  low,  St  can 
be  as  large  as  0.2,  but  in  most  of  our  experiments  gravitational 
effects  are  negligible.  (Note:  Centrifugal  effects  are  always 
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negligible).  Obviously,  a  less  restrictive  model  of  the  wiping 
process  is  required,  but  since  an  analytical  solution  is  very 
unlikely,  it  is  necessary  to  consider  the  use  of  numerical 
techniques. 

We  are  presently  attempting  to  use  the  Finite  Element  Method 
(FEM)  to  solve  the  blade  wiping  problem  for  newtonian  fluids 
without  the  simplifying  assumptions  mentioned  above.  Surface 
tension,  which  is  neglected  by  the  use  of  the  Swif t-steiber 
boundary  conditions  in  the  lubrication  theory  model,  can  be 
accounted  for  by  treating  the  region  downstream  of  the  blade  as 
a  free  surface  problem.  The  effects  of  gravity,  inertia,  and 
curvature  are  relatively  easy  to  incorporate  into  the  numerical 
scheme  as  well,  so  that  the  rigid  blade  wiping  problem  can  be 
completely  solved  for  newtonian  liquids. 

Even  with  sophisticated  numerical  techniques,  it  is  unlikely 
that  the  wiping  problem  can  be  solved  in  general  for  viscoelastic 
liquids.  With  newtonian  liquids,  the  formulation  of  the  problem 
is  straightforward  because  of  the  Bimple  and  predictable  manner 
in  which  they  respond  to  imposed  deformations.  Viscoelastic 
liquids,  on  the  other  hand,  exhibit  extremely  complex  behavior 
which  is  dependent  on  both  the  magnitude  and  duration  of 
deformation  imposed.  Unlike  the  newtonian  case,  a  constitutive 
equation  does  not  exist  that  is  capable  of  describing  visco¬ 
elastic  fluid  behavior  in  a  general  flow  field.  Nonlinear 
viscous  behavior  alone,  is  relatively  easy  to  model,  the  power 
law  fluid  being  a  simple  example.  The  difficulties  arise  because 
the  viscous  effects  are  coupled  with  transient  and  deformation 
history  dependent  behavior  associated  with  the  "elasticity"  of 
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the  fluid.  Today,  flow  problems  are  solved  using  constitutive 
equations  chosen  according  to  a  compromise  between  algebraic 
tractability  and  physical  reality.  Such  equations  are  highly 
dependent  on  the  nature  of  the  flow  field  and  are  limited  to 
special  classes  of  problems.  The  wiping  problem  is  a  formidible 
problem  for  a  constitutive  equation  to  model  because  fluid 
particles  experience  rapidly  changing  deformations  in  time,  which 
can  be  of  considerable  importance  for  viscoelastic  fluids 
(Shirodkar,  1982).  Elastic  liquids  can  exhibit  large  transient 
shear  stresses  as  well  as  extensional  stresses  when  undergoing 
the  rapid  changes  in  deformation  rate  that  are  present  because  of 
the  converging  nature  of  the  flow.  The  reulting  kinematics  for 
viscoelastic  liquids  are  often  drastically  different  from  those 
present  for  newtonian  liquids  (Doremus,  1983,  Mensah,  1984, 
Metzner,  (1969). 

Key  to  solving  the  problem  of  modeling  viscoelastic  fluid 
behavior  is  determing  which  material  properties  are  most 
responsible  for  the  behavior  of  a  given  fluid  in  a  particular 
flow  field.  We  have  mentioned  the  recoverable  shear,  which  is  a 
steady  viscometric  material  property  providing  a  measure  of  the 
relative  importance  of  elastic  and  viscous  effects  in  steady 
shear  flow.  An  important  question  is,  can  material  properties 
measured  in  steady  viscometric  flows  provide  insight  into  how  an 
elastic  liquid  will  behave  in  a  transient  flow  field?  If  not, 
what  are  the  pertinent  transient  fluid  properties  to  measure  and 
how  can  they  be  incorporated  into  useable  and  realistic 
constitutive  equations? 


Given  the  difficulty  involved  in  choosing  an  appropriate 
constitutive  equation  to  model  the  wiping  process,  it  is 
fortunate  that  we  are  able  to  use  experimental  data  to  evaluate 
proposed  constitutive  equations.  Finding  a  constitutive  equation 
to  best  fit  the  data  is  equivalent  to  determining  which  material 
properties  of  the  fluid  dominate  the  flow  field.  As  a  simple 
example,  if  the  power-law  model  had  fit  the  data  nicely,  we  would 
conclude  that  the  shear  thinning  nature  of  the  fluid  is  most 
responsible  for  the  resulting  behavior.  In  fact,  such  is  the  case 
for  the  CMC  solution  with  the  lowest  recoverable  shear.  There  is 
an  increase  in  X  over  newtonian  values  as  predicted  by  the  model, 
although  a  little  more  dramatic  than  expected.  However,  as  the 
elasticity  of  the  fluid  is  increased  by  the  addition  of  poly¬ 
acrylamide,  as  reflected  in  an  increased  recoverable  shear,  the 
trend  in  X  is  reversed.  The  shear  dependent  viscosity  is 
obscured  by  another  fluid  property  associated  with  the  elasticity 
of  the  fluid.  The  data  provide  experimental  evidence  of  this 
fact,  and  is  a  good  test  for  evaluating  constitutive  equations  to 
be  used  to  model  the  process.  In  this  way  we  can  obtain  a  better 
understanding  of  the  important  features  of  the  wiping  process. 

This  study  provides  a  good  foundation  for  understanding  the 
rigid  blade  wiping  process,  but  equally  important,  it  provides 
some  insight  into  a  direction  for  future  work.  A  Finite  Element 
solution  to  the  newtonian  problem  is  required  to  understand  the 
effects  of  inertia  and  surface  tension  and  explore  thin  wipers 
and  large  wiper  angle  geometries.  With  the  newtonian  problem 
well  understood,  experimental  data  for  well  characterized  visco¬ 
elastic  fluids  can  be  compared  to  newtonian  results  and  provide 
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information  regarding  more  quantitative  aspects  of  non-newtonian 
and  elastic  effects.  An  obvious  extension  of  the  rigid  blade 
wiping  problem  is  to  the  study  of  wiping  with  flexible  blades. 
This  problem  is  much  more  complicated  than  the  rigid  wiper  case 
because  the  geometry  is  unknown  a-priori  and  depends  upon  the 
complex  interaction  of  hydrodynamic  forces  and  wiper  properties 
(Saita,  1984).  Fortunately,  the  solution  to  the  rigid  wiper 
problem  provides  the  logical  first  step  in  understanding  the 
influence  of  geometry  on  the  hydrodynamic  forces  that  develop 
under  a  wiper  and  is  a  valuable  tool  for  further  study. 


tr  I 


F/W  -  force  per  unit  width  on  blade 
h  -  residual  fluid  thickness 

-  characteristic  length  in  the  y-direction 
H  -  blade  height  at  x«0 

H0  -  minimum  clearance  between  the  blade  and  the  roll 
K  -  parameter  in  power-law  fluid  model 

L  -  blade  length 

L*  -  entrance  length  (Appendix  A) 

Lx  -  characteristic  length  in  x-direction 

n  -  power-law  index 

N1  -  primary  normal  stress  coefficient  (page  30) 

P  -  pressure 

Q  -  flow  rate 

R  -  roll  radius 

Re  -  Reynolds  number  (page  5) 

U  -  linear  roll  speed 

u  -  velocity  in  x-direction 

u (x)  -  u-velocity  on  roll  surface 

v  -  velocity  in  y-direction 

v(x)  -  v-velocity  on  roll  surface 

x  -  space  coordinate  in  the  primary  flow  direction 

X3  -  x-coordinate  at  the  point  of  minimum  blade  clearance 

XI  -  x-coordinate  at  the  upstream  edge  of  the  blade 

X2  -  x-coordinate  at  the  downstream  edge  of  blade 

y  -  space  coordinate  transverse  to  primary  flow  direction 


Variables: (page  2) 


(*'<x) 
p  (5) 


S  (x) 


A  - 


7 
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5 

TT 

a. 


LA 

f 

cr'(x) 


a  (?) 

, « 


CT 

r 


x.. 

*j 


5 

<AJ 


-  geometric  parameter  (page  5) 

-  function  describing  the  blade  surface  (linear) 

-  dimensionless  blade  surface 

-  boundary  layer  thickness  (Appendix  A) 
rate  of  deformation  tensor 

-  dimensionless  space  coordinate  -  y/h  (page  4) 

-  dimensionless  residual  fluid  thickness  -  h/h  (page  6) 

-  viscosity 

-  dimensionless  space  coordinate  -  x/LA  (page  4) 

-  dimensionless  pressure  (page  4) 

"  second  invariant  of  A 

-  density 

-  function  describing  the  roll  surface  (quadratic) 

-  dimensionless  roll  surface 

-  surface  tension 

-  stress  tensor 

-  i,j  component  of  the  stress  tensor 

-  dimensionless  u-velocity  (page  4) 

-  dimensionless  u(x)  (page  11) 

-  dimensionless  blade  loading  (page  13) 

-  dimensionless  v-velocity  (page  4) 

-  dimensionless  v(x)  (page  11) 

-  blade  angle 
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FIGURE  It  Conf igurat ion  of  a  simple 
Blada-Over-Rol 1  Wiper 


FIGURE  2j  Uafinition  sketch  of  the  Geometry 
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FIGURE  3i  Effact  of  blada  angle  on  Residual 

Fluid  Thicknoss  C  ^  5,  Solid  linos 
ara  Lubrication  Theory  for  a  Newtonian 
liquid.  Symbols  are  from  Hsu  <1984) 
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FIGURE  4i  Effect  of  blade  displacement  on  Residual 
Fluid  Thickness  <  A  ).  Solid  lines 
are  Lubrication  Theory  for  a  Newtonian 
liquid.  Symbols  are  from  Hsu  <1984). 

O  -  Zeroth  order  solution 
®  -  First  order  curvature  included 
^  -  Finite  Element  Solution 
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FIGURE  5«  Effect  of  blade  angle  on  Blade  Loading 
Solid  line#  ore  Lubrication  Theory 
for  a  Newtonian  liquid.  Symbols  are 
from  Hsu  (19845. 

0  -  Zeroth  order  solution 

I  -  First  order  curvature  included 
A,  -  Finite  Element  Solution 
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FIGURE  6<  Effect  of  blada  displacement  on  Blade 
Loading.  Solid  lines  ore  Lubrication 
Theory  for  a  Newtonian  liquid. 

Symbols  are  from  Hsu  (1904). 
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Effect  of  blade  length  on  Blade  Loading 
Solid  lines  are  Lubrication  Theory 
for  a  Newtonian  liquid. 
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FIGURE  9ai  Comparison  of  the  Pressure  Distribution 

under  ths  blads  obtained  with  the  numerical 
Power- law  Lubrication  Theory  algorithm 
(n-1)  to  the  analytical  Newtonian 
Lubrication  Theory  solution. 
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FIGURE  9b«  Comparison  of  the  Residual  Fluid 

Thickness  obtained  with  the  numerical 
Powei — law  Lubrication  Theory  algorithm 
Cn«l)  to  the  analytical  Newtonian 
Lubrication  Theory  solution. 
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FIGURE  IQs  Effect  of  powei — law  index  on  Residual 
Fluid  Thickness  with  blade  angle  as  a 
parameter.  (  Key  in  Table  1  > 
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FIGURE  lli  Effect  of  powei — law  index  on  Residual 
Fluid  Thickness  with  blade  angle  us  a 
parameter.  C  Key  in  Table  1  > 
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FIGURE  13*  Effect  of  Power-law  index  on  Blade  Loading 
with  blade  angle  and  blade  height  as 
parameters.  (  Key  in  Table  1  ) 
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FIGURE  I4i  Effect  of  Power-law  index  on  Residual 

Fluid  Thickness  with  blade  displacement 
as  a  parameter.  C  Key  in  Table  1  > 
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FIGURE  35«  Effoct  of  Power-law  index  on  Reeidual 

Fluid  Thickneee  with  blade  displacement 
□e  a  parameter.  C  Key  in  Table  1  > 
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FIGURE  17i  Effect  of  Powar-law  index  on  Blade  Loading 
with  blade  displacement  and  blade  height 
as  parameters.  C  Key  in  Table  1  ) 


FIGURE  19«  Blade  Geometry  -  Length 
and  Angle  <  w  )  are  adjustable. 

CNotei  L  Is  the  land  length] 
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FIGURE  20i  Rasidual  Fluid  Thickness  vs.  Blade  Height 
for  Glycerine  <o)  and  Kara  Syrup  (x) 
compared  tc  Lubrication  Theory  Cline) 


Cl/1000  in) 
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FIGURE  21s  Newtonian  data  and  Lubrication  Theory 
showing  tha  offset  of  blade  longth. 
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FIGURE  22i  Newtonian  data  and  Lubrication  Theory 
showing  the  effect  of  blade  angle 
using  a  1/2  in.  blade. 
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FIGURE  23j  Newtonian  data  and  Lubrication  Theory 
showing  the  effect  of  blade  angle 
using  a  1/4  in.  blade. 
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Effect  of  Slade  Displacement  <X2> 
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FIGURE  25i  Newtonian  data  and  Lubrication  Theory 

showing  the  effect  of  blade  displacement 
using  a  1/2  in.  blade  at  0  degrees. 
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Effect  of  Shear-thinning  and  Elastic  Behavior 


R-2.  75  in 
L-0.  25  in 
X2-0 

10  deg  angle 

Glycerine 
1.25%  CMC 

1.25%  CMC  and  0.02%  PAM 
1.25%  CMC  and  0. 03%  PAM 


Solid  lines  -  Lubrication  Theory 
Uppen  Power- law  with  n-0.  8 
Loweri  Newtonian 


Ho  Cl/1000  in) 


FIGURE  28s  Data  for  Newtonian  and  Viscoelastic 
liquids  compared  to  Newtonian  and 
Powei — law  Lubrication  Theory. 
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FIGURE  32»  Data  for  Newtonian  and  Viscoelastic 
liquids  compared  to  Newtonian  and 
Powei — law  Lubrication  Theory. 
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FIGURE  33i  Data  for  Newtonian  and  Viscoelastic 
liquids  compared  to  Newtonian  and 
Powei — law  Lubrication  Theory. 
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APPENDIX  A 


MOMENTUM  TRANSFER  ARGUMENTS  FOR  AN  INCREASE  IN  RESIDUAL  FLUID 
THICKNESS  WITH  DECREASING  BLADE  LENGTH. 

For  the  purpose  of  simplifying  the  present  discussion,  we 
consider  the  simplified  geometry  with  zero  curvature  (R-#>  )  and 
zero  wiper  blade  angle  shown  in  Figure  A-l.  Similar  arguments 
apply  to  more  complicated  geometries. 

Fluid  entrained  by  the  moving  roll  in  our  experimental 
apparatus  approaches  the  wiper  at  a  constant  speed,  equal  to  the 
roll  speed.  However,  because  of  the  no-slip  condition  at  solid 
surfaces,  the  velocity  of  fluid  making  contact  with  the  leading 
edge  of  the  wiper  becomes  zero,  which  initiates  the  formation  of 
a  boundary  layer  as  shown  in  Figure  A-2,  The  thickness  of  the 
boundary  layer  ,  J(x),  is  the  distance  over  which  viscous  shear 
stresses,  originating  at  the  wiper  surface,  transfer  momentum 
between  the  wiper  and  the  fluid.  Fluid  outside  of  the  boundary 
layer  remains  unaffected  by  the  presence  of  the  wiper.  The 
boundary  layer  grows  along  the  wiper,  downsteam  of  the  leading 
edge,  and  if  the  wiper  is  long  enough,  eventually  fills  the 
region  between  the  wiper  and  the  moving  surface.  For  two 
parallel  surfaces  ,  one  moving  and  the  other  stationary,  the 
velocity  profile  is  linear  when  it  becomes  "fully  developed",  and 
the  distance  required  for  this  to  occur  is  called  the  "entrance 
length  (Lft)".  Simple  lubrication  theory  assumes  that  the 
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entrance  length  is  zero  and  therefore,  cannot  account  for 
"entrance  effects",  should  they  occur  due  to  the  developing 
boundary  layer.  For  large  wipers,  when  L/L&  >>  1,  entrance 
effects  are  negligible,  but  as  the  wiper  length  decreases,  the 
entrance  region  can  become  a  significant  fraction  of  the  wiper 
length.  We  would  like  to  have  some  estimate  of  how  small  a  wiper 
must  be  in  order  for  entrance  effects  to  be  significant  and  how 
these  effects  can  influence  the  residual  fluid  thickness. 

In  an  attempt  to  investigate  entrance  effects  we  consider 
the  simple  geometry  shown  in  Figue  A-l.  From  boundary  layer 
theory  (Bird,  1960  -  pg.145)  we  know  that  the  thickness  of  a 
boundary  layer  developing  on  the  wiper  can  be  approximated  by 

<5(  =  H,  J  “  (a-D 


with  the  velocity  distribution 


5 ;  *  (&)■*(  fa*) 

j  i 


(A-2) 


Assuming,  that  outside  of  the  boundary  layer  fluid  moves  with  the 
same  velocity  as  the  moving  surface  (u»U) ,  the  flow  rate,  Q, 
between  the  wiper  and  the  moving  surface  can  be  written 


(A- 3) 


where  &(x)  <  H  .  If  the  boundary  layer  thickness  in  the  above 
0 

expression  is  evaluated  at  the  downstream  edge  of  the  wiper  (i.e. 


79 


x=L) ,  then  the  residual  fluid  thickness  is  given  by 


(A- 4) 


Using  the  following  dimensionless  variables  and  substituting 


Hd  ^  P 


Equations  A-l  and  A-4  into  A-3,  we  obtain  an  expression  for  the 
residual  fluid  thickness/  which  is  only  valid  as  long  as  the 
boundary  layer  does  not  extend  to  the  moving  surface  (i.e. 
J(L)  *  H0). 


(A-5) 


The  assumption  that  <5  (L)  v<  H  constrains  the  blade  length  to 

.  «,  I40^ 

L  '  (M.bM'l1  (A‘6) 

The  actual  wiper  lengths  for  which  entrance  effects  are 
significant  is  actually  larger  than  this,  because  once  the 
boundary  layer  reaches  the  moving  surface,  there  is  an  additional 
transition  region  along  which  the  velocity  transforms  into  its 
fully  developed  linear  profile. 

Figure  A-3  shows  the  effect  of  blade  length,  at  a  given 
Reynolds  number,  on  the  residual  fluid  thickness  as  predicted  by 


0 


8 


Equation  A-5.  According  to  the  boundary  layer  arguments,  the 
residual  fluid  thickness  increases  as  the  blade  length  decreases, 
when  entrance  effects  are  significant.  The  dotted  line  is  an 
extension  of  Equation  A-5  into  the  transition  region  of  blade 
length,  where  boundary  layer  arguments  are  certainly  invalid. 
Although  we  can  not  easily  estimate  the  length  of  the  transition 
region,  it  is  obvious  that  simple  momentum  transfer  arguments 
suggest  that  the  residual  fluid  thickness  increases  when  entrance 
effects  are  significant. 


81 


Figure  A-1j  Simple  wiper  geometry  with  zero 
curvature  and  zero  wiper  angle. 


Figure  A-2:  Developing  boundary  layer  at 
the  leading  edge  of  a  wiper. 
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Figure  A-3i 


Dependence  of  residual  fluid 
thickness  on  blade  length  when 
entrance  effects  due  to  developing 
boundary  layer  are  significant. 
Specific  for  geometry  in  Figure  A-l 
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